
import delimited "Additional_controls_df.csv", clear 

encode iso3, generate(ID)
xtset ID year

g ColdWar=0
replace ColdWar=1 if year<=1991
label var ColdWar "Cold War"

g l_msc = l.months_since_coup
label var l_msc "Months since coup"
g l_msc_sq = l_msc*l_msc
label var l_msc_sq "Months since coup (sq)"
g l_msc_cu = l_msc*l_msc*l_msc
label var l_msc_cu "Months since coup (cu)"

g l_ncp5 = l.n_coup_5y
label var l_ncp5 "N. coups 5 years (cu)"

g l_ncp5_sq = l_ncp5*l_ncp5
label var l_ncp5_sq "N. coups 5 years (cu)"

g l_ncp5_cu = l_ncp5*l_ncp5*l_ncp5
label var l_ncp5_cu "N. coups 5 years (cu)"


label var colony_esp "Colony Esp"
label var colony_gbr "Colony GBR"
label var colony_fra "Colony FRA"
label var colony_prt "Colony PRT"
label var colony_oeu "Colony OEU"


probit pt_attempt   c.l.months_since_coup##c.l.months_since_coup##c.l.months_since_coup ///
 c.l.n_coup_5y##c.l.n_coup_5y##c.l.n_coup_5y ///
 i.colony_esp i.colony_gbr i.colony_fra i.colony_prt i.colony_oeu i.year ColdWar ,   vce(cluster ID) 

lroc, nograph

predict CoupRiskPr, p 


probit pt_attempt   l_msc l_msc_sq l_msc_cu ///
 l_ncp5 l_ncp5_sq l_ncp5_cu ///
 i.colony_esp i.colony_gbr i.colony_fra i.colony_prt i.colony_oeu i.year ColdWar ,   vce(cluster ID) 

outreg2 using Tab1.doc, ctitle(Model 1, Coup, Probit) ///
sortvar ///
(  l_msc l_msc_sq l_msc_cu l_ncp5 l_ncp5_sq l_ncp5_cu colony_esp colony_gbr colony_fra colony_prt colony_oeu  ColdWar  ///
) ///
keep ///
(  l_msc l_msc_sq l_msc_cu l_ncp5 l_ncp5_sq l_ncp5_cu colony_esp colony_gbr colony_fra colony_prt colony_oeu  ColdWar  ///
) ///
tex dec(3) pdec(3) addtext(Year FE, YES, Cluster SE, YES) label  replace  



 
save "Additional_controls_crp_df.dta", replace


import delimited "morning_after_dataset_countrylevel_additional_controls.csv", clear 


g COUNTRY_ID= country_id
g YEAR_DATA= year_data


********************************************************************************
***Appendix A. Replicates Figure 2 in the paper***
********************************************************************************
estimates clear 

***Table A1*** coupr risk as measure
*Variation to Figure 2 in the article (left panel, Country FE)*
xtset COUNTRY_ID YEAR_DATA
*Base model
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store pl_c_b
*Main model
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ma_c_b
*Electoral controls
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store el_c_b
*Economic controls
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.oil_valuepop_2014_ln L.growth_pwt i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ec_c_b
*Instability controls
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store in_c_b
*Purges MAs
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ym_c_b
*All controls
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a L.oil_valuepop_2014_ln L.growth_pwt L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ac_c_b

***Table A2*** coupr risk as measure
*Variation to Figure 2 in the article (left panel, Leader FE)*
xtset leader_cab_id YEAR_DATA
*Base model
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store pl_l_b
*Main model
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ma_l_b
*Electoral controls
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store el_l_b
*Economic controls
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.oil_valuepop_2014_ln L.growth_pwt i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ec_l_b
*Instability controls
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store in_l_b
*Purges MAs
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ym_l_b
*All controls
xtreg replacement_rateadj_minister coupriskpr L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a L.oil_valuepop_2014_ln L.growth_pwt L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ac_l_b


********************************************************************************
***Figure 2 graphic***
********************************************************************************

*Requires the Coefplot packacge

ssc install coefplot, replace

graph set window fontface "Times New Roman" 
graph set window fontfacemono "Times New Roman"
graph set window fontfacesans "Times New Roman"
graph set window fontfaceserif "Times New Roman"


*Replicate graphic
coefplot (pl_c_b, offset(0.45) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(pl_l_b, offset(0.43) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ma_c_b, offset(0.30) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ma_l_b, offset(0.28) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(el_c_b, offset(0.15) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(el_l_b, offset(0.13) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ec_c_b, offset(0.0) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ec_l_b, offset(-0.02) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(in_c_b, offset(-0.15) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(in_l_b, offset(-0.17) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ym_c_b, offset(-0.30) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ym_l_b, offset(-0.32) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ac_c_b, offset(-0.45) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ac_l_b, offset(-0.47) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))), ///
keep(coupriskpr) yscale(noline alt) grid(glcolor(gs6) glwidth(thin)) ///
xline(0, lcolor(black) lpattern(shortdash)) ytitle("", size(small)) levels(95) title("{bf:Failed coup attempts - coup risk as IV}", color(black) margin(medsmall) size(medium)) ///
graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) coeflabels(, nolabel notick) ///
legend(on region(lstyle(none) lcolor(white)) order(6 4) label(6 "Country FE") label(4 "Leader FE") size(*.75) symxsize(*.3) symysize(*0.8) colgap(2) keygap(1) forcesize ) ///
xsize(12) ysize(18) xlabel(-0.05(0.5)1.55, labsize(vsmall)) ylabel(, labsize(vsmall)) baselevels msymbol(o) ///
text(0.55 1.45 "Base models", size(small)) ///
text(0.70 1.45 "Main models", size(small)) ///
text(0.85 1.45 "Electoral controls", size(small)) ///
text(1.00 1.45 "Economic controls", size(small)) ///
text(1.15 1.45 "Instability controls", size(small)) ///
text(1.30 1.45 "Purges MAs", size(small)) ///
text(1.45 1.45 "All controls", size(small)) name(gr1a) 
graph save gr1a, replace


****************************** original models *********************************
********************************************************************************
***Appendix A. Replicates Figure 2 in the paper***
********************************************************************************
estimates clear 

***Table A1***
*Equivalent to Figure 2 in the article (left panel, Country FE)*
xtset COUNTRY_ID YEAR_DATA
*Base model
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store pl_c
*Main model
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ma_c
*Electoral controls
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store el_c
*Economic controls
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.oil_valuepop_2014_ln L.growth_pwt i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ec_c
*Instability controls
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store in_c
*Purges MAs
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ym_c
*All controls
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a L.oil_valuepop_2014_ln L.growth_pwt L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ac_c

***Table A2***
*Equivalent to Figure 2 in the article (left panel, Leader FE)*
xtset leader_cab_id YEAR_DATA
*Base model
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store pl_l
*Main model
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ma_l
*Electoral controls
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store el_l
*Economic controls
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.oil_valuepop_2014_ln L.growth_pwt i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ec_l
*Instability controls
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store in_l
*Purges MAs
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ym_l
*All controls
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a L.oil_valuepop_2014_ln L.growth_pwt L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ac_l


*Replicate graphic
coefplot (pl_c, offset(0.45) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(pl_l, offset(0.43) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ma_c, offset(0.30) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ma_l, offset(0.28) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(el_c, offset(0.15) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(el_l, offset(0.13) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ec_c, offset(0.0) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ec_l, offset(-0.02) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(in_c, offset(-0.15) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(in_l, offset(-0.17) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ym_c, offset(-0.30) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ym_l, offset(-0.32) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ac_c, offset(-0.45) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ac_l, offset(-0.47) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))), ///
keep(coupattempt_whogov) yscale(noline alt) grid(glcolor(gs6) glwidth(thin)) ///
xline(0, lcolor(black) lpattern(shortdash)) ytitle("", size(small)) levels(95) title("{bf:Failed coup attempts - original IV}", color(black) margin(medsmall) size(medium)) ///
graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) coeflabels(, nolabel notick) ///
legend(on region(lstyle(none) lcolor(white)) order(6 4) label(6 "Country FE") label(4 "Leader FE") size(*.75) symxsize(*.3) symysize(*0.8) colgap(2) keygap(1) forcesize ) ///
xsize(12) ysize(18) xlabel(-0.05(0.05)0.35, labsize(vsmall)) ylabel(, labsize(vsmall)) baselevels msymbol(o) ///
text(0.55 0.27 "Base models", size(small)) ///
text(0.70 0.27 "Main models", size(small)) ///
text(0.85 0.27 "Electoral controls", size(small)) ///
text(1.00 0.27 "Economic controls", size(small)) ///
text(1.15 0.27 "Instability controls", size(small)) ///
text(1.30 0.27 "Purges MAs", size(small)) ///
text(1.45 0.27 "All controls", size(small)) name(gr1b) 
graph save gr1b, replace


graph combine gr1a gr1b, xsize(20) ysize(14) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) rows(1) cols(2) imargin(medium) name(gr1) 


***Table A1*** coup proofing as additional control
*Equivalent to Figure 2 in the article (left panel, Country FE)*
xtset COUNTRY_ID YEAR_DATA
*Base model + de bruing coup proofing data
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln l.cbcount l.ha_cbcount l.pcount l.effectivenumber i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store pl_c_c
*Main model + de bruing coup proofing data
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c l.cbcount l.ha_cbcount l.pcount l.effectivenumber i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ma_c_c
*All controls + de bruing coup proofing data
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a L.oil_valuepop_2014_ln L.growth_pwt L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 l.cbcount l.ha_cbcount l.pcount l.effectivenumber i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ac_c_c

***Table A1*** coup proofing as additional control
*Equivalent to Figure 2 in the article (left panel, Leader FE)*
xtset leader_cab_id YEAR_DATA
*Base model + de bruing coup proofing data
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln l.cbcount l.ha_cbcount l.pcount l.effectivenumber i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store pl_l_c
*Main model + de bruing coup proofing data
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c l.cbcount l.ha_cbcount l.pcount l.effectivenumber i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ma_l_c
*All controls + de bruing coup proofing data
xtreg replacement_rateadj_minister coupattempt_whogov L.gdp_cap_pwt_ln L.pop_pwt_ln L.military_c L.monarchy_c L.party_c L.elec_all L.t1_e_a L.oil_valuepop_2014_ln L.growth_pwt L.onset1 L.nonvio_camp_2 L.mid_onset L.war_onset L.domestic3_10 replacement_rateadj_minister_300 replacement_rateadj_minister_200 replacement_rateadj_minister_100 l.cbcount l.ha_cbcount l.pcount l.effectivenumber i.YEAR_DATA if democracy1==0 & independent==1, fe cluster(COUNTRY_ID)
estimates store ac_l_c

coefplot (pl_c_c, offset(0.45) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(pl_l_c, offset(0.43) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ma_c_c, offset(0.05) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ma_l_c, offset(0.08) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))) ///
(ac_c_c, offset(-0.35) color(black) ciopts(lcolor(black black) lwidth(thin medthick))) ///
(ac_l_c, offset(-0.38) color(gs10) ciopts(lcolor(gs10 gs10) lwidth(thin medthick))), ///
keep(coupattempt_whogov) yscale(noline alt) grid(glcolor(gs6) glwidth(thin)) ///
xline(0, lcolor(black) lpattern(shortdash)) ytitle("", size(small)) levels(95) title("{bf:Failed coup attempts - CP controls}", color(black) margin(medsmall) size(medium)) ///
graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) coeflabels(, nolabel notick) ///
legend(on region(lstyle(none) lcolor(white)) order(6 4) label(6 "Country FE") label(4 "Leader FE") size(*.75) symxsize(*.3) symysize(*0.8) colgap(2) keygap(1) forcesize ) ///
xsize(6) ysize(4) xlabel(-0.05(0.05)0.35, labsize(vsmall)) ylabel(, labsize(vsmall)) baselevels msymbol(o) ///
text(0.55 0.27 "Base models", size(small)) ///
text(0.95 0.27 "Main models", size(small)) ///
text(1.35 0.27 "All controls", size(small)) name(gr2a) 





